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Abstract. Adopting a purely two dimensional relativistic equation for graphene's 

carriers contradicts the Heisenberg uncertainty principle since it requires setting off- 
the-surface coordinate of a three-dimensional wavefunction to zero. Here we present 
a theoretical framework for describing graphene's massless relativistic carriers in 
accordance with this most fundamental of all quantum principles. A gradual confining 
procedure is used to restrict the dynamics onto a surface and normal to the surface 
parts and in the process the embedding of this surface into the three dimensional 
world is accounted for. As a result an invariant geometric potential arises in the 
surface part which scales linearly with the Mean curvature and shifts the Fermi energy 
of the material proportional to bending. Strain induced modification of the electronic 
properties or "straintronics" is clearly an important field of study in graphene. This 
opens a venue to producing electronic devices, MEMS and NEMS where the electronic 
properties are controlled by geometric means and no additional alteration of graphene 
is necessary. The appearance of this geometric potential also provides us with clues 
as to how quantum dynamics looks like in the curved space-time of general relativity. 
In this context, we explore a two-dimensional cross-section of the wormhole geometry 
realized with graphene as a solid state thought experiment. 
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1. Introduction 

The two-dimensional allotrope of carbon, namely graphene, has emerged as a very 
promising electronic material. In monolayer graphene electrons are described by the 
Dirac equation [T] . Here we analyze the geometrically induced chemical potential due to 
confinement (i.e. the gauge field) in the Dirac equation on a two dimensional sheet since 
graphene is a two dimensional sheet. Height fluctuations of a graphene layer (carriers 
are described by a relativistic equation) on a rough substrate generate a nonvanishing 
geometric potential. In this case, the shape of the graphene sheet is determined by a 
competition between the interaction of the layer with the rough substrate Hs^b (which 
tends to impose a preferred height), the elastic properties of the layer Heiast and the 
electronic energy ifgi, which is also a function of geometry: 

H = Hsub + Heiast + Hei. 

The conflning procedure employed here for exploring the electronic properties of 
graphene goes beyond the standard analysis of the Dirac equation in curved geometries 
^ . Here we will employ an alternative approach which is consistent with the Heisenberg 
principle. On one hand the starting point is the realization that graphene as a one- 
atom-thick membrane has carriers conflned in a two dimensional space trapped in three 
dimensions. Furthermore the electrons are described by a massless relativistic equation. 
On the other hand the wavefunction of a quantum particle is always three dimensional 
due to the Heisenberg uncertainty principle which forbids setting any of the coordinates 
to zero [this would lead to indefiniteness in the momentum according to Ap > h/ {2 Ax)]. 
In this way any two dimensional quantum motion would have an evanescent off- 
surface component of the wavefunction which can probe the two dimensional surface 
for curvature [21 S]. If the graphene sheet is curved, either intrinsically or extrinsically, 
then the carriers will at least be able to "feel" this curvature in the form of some mass 
or effective gauge field. 

In fact, curvature of the sheets builds strain and from a microscopic point of view 
strain has been shown to modify the electronic structure [5]. This is a confirmation of 
the line of thinking presented in this paper. When interatomic distances are modified 
due to strain in the underlying lattice caused by curvature, the periodicity is disrupted 
and the conditions for the application of Bloch theorem do not hold. Therefore, standard 
methods of solution fail to grasp the complexity of the rippling property of graphene. 
Effectively this means that one has to use a constraining procedure which starts from 
three dimensions and gradually confines the quantum dynamics onto the surface. Note, 
for the Schrodinger equation the curvature induces an effective (da Costa) potential [1] 
through a confining procedure reducing the dimensionality of the quantum system. One 
question stands out: what is the corresponding potential for the relativistic case? We 
have answered this question and illustrated it through graphene. In short, the effective 
potential has a linear dependence on the Mean curvature (M) as opposed to — K, 
where K is the Gaussian curvature, in the usual (or nonrelativistic) case. 
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This effect is studied by expanding around the two dimensional space of the 
graphene sheet for vanishingly small excursions in the third direction. Such a procedure 
is well known for the Schrodinger equation [3] and we have carried it out in the case 
of a relativistic equation. The idea behind this confinement is that an "external force" 
gradually compresses the quantum dynamics onto a surface. In the following section we 
discuss where this force appears from in the case of graphene. When the excursion in 
the normal to the surface direction becomes small enough one can take the limit (the 
normal to the surface coordinate goes to zero) and split the original 3+1 dimensional 
relativistic equation into 2+1 and 1+1 dimensional equations. The 2+1 dimensional 
equation encodes the quantum dynamics onto the surface and the 1+1 dimensional 
equation describes the behavior of the "evanescent component", that is the normal to 
the surface component of the wavefunction. In this way the Heisenberg principle is not 
violated since in the process described above no setting-to-zero of any coordinate of a 
three dimensional wavefunction takes place. 

In reducing the Dirac equation from 3+1 dimensions to 2+1 and 1+1 dimensional 
equations we use the fundamental 2 representation of the Clifford algebra |6], which is 
appropriate for the isospin of graphene. The results of this confining procedure (Sec. II) 
and its specific applications (Sec. Ill) are also relevant for other fields of physics where 
phenomena are described within the framework of Dirac theory: chiral spin liquids, 
quantum Hall systems, anyons. 

Experimental confirmations, as to what is the correct relativistic description of 
quantum motion on a curved manifold, are important. Graphene presents with an 
experimental opportunity to answer this question with potential applications to the 
theory of gravity and its connection to quantum mechanics. This aspect also motivates 
us to consider the two-dimensional cross-sections of a wormhole geometry (Sec. IV). 
Our main findings are summarized in Sec. V. 

2. Confinement 

Now let us consider constrained nonrelativistic dynamics described by the Schrodinger 
equation with the extra chemical potential due to confinement^ 



where Ki are the principal curvatures of the surface, h is the Planck's constant and m 
is the effective mass. Here M is the Mean curvature and K is the Gauss curvaure. 
Note, by virtue of its derivation, this potential stems from the kinetic part of the 
constrained Schrodinger equation, that is in the case of the massless Dirac equation 
geometric potentials stemming from the kinetic part will not vanish. 

As is well known, the Dirac equation is a "square root" of the Schrodinger equation, 
meaning that whatever the constrained Dirac equation is, we would expect the chemical 
potential to scale linearly with the curvature since it scales as the square of the curvature 
in the case of the Schrodinger equation. 




(1) 
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Figure 1. Bending the surface of grapliene by a radius R and its effect on the Pz 
orbitals. 



Consider bending the graphene sheet by a radius R. This will decrease the distance 
between the orbitals from / = R6 to d = 2i?sin (|), 

4. --2^«l- P) 
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The decrease in the distance between the orbitals increases the overlap between the two 
lobes of the p^ orbital, see Fig. [T| The overlap scales as the square of the curvature 
\/R but this does not necessarily point to an effective chemical potential proportional 
to the square of the curvature in the Dirac equation for the carriers[71 [1]. 

Next, let us proceed with the derivation. Here we choose to use gauge covariant 
Dirac description for the massless carriers in graphene, both electrons and holes: 

-2nt;ir^7'^D, + m(g3)4 + ^±fe)] ^(t,gi) = 0, (3) 

where (gi, 52) are the generalized coordinates of the surface and measures the deviation 
from the surface in the embedding 3 dimensional space. Here Da stands for a 3+1 
dimensional gauge covariant derivative, 7" are curved gamma matrices 7** = 7"(q'i, q2i Qs) 
related to the usual constant matrices 7^ by 

T^l^i, 92, gs) = e^^iqi, q2, q3)Y, (4) 

where is the tetrad constructed from the metric tensor. Here m{q3) is an effective 
mass term with the following property 

."(*) = { ° * = ° (5) 

and V±{q3) is the constraining potential forcing the system onto the surface 

Here a possible model of the confining potential is the charged plate potential 

V±iQ3) ~ Igsl (7) 

since when the electron deviates from the two dimensional world of the sheet it disrupts 
the charge balance and a Coulomb force emerges. 
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The gauge covariant derivative is defined through its action on a vector field as 
follows 

Do = du = d^v' + Tly - %-A,v\ (8) 

c 

where Aj are the components of a vector potential and are the Christoffel symbols 
of the second kind = ^G^"^ {Gmk,i + Gmi,k — Gki,m) ■ Here Gmk is the metric tensor 
of the embedding space. This equation is defined in the entire embedding space 

Our goal now is to constrain ^ onto a surface x{qi,q2) and separate in two 
parts: one that describes the relativistic dynamics on the surface and one that 
describes the dynamics in a direction q^ orthogonal to the surface. We will treat 
the off surface dynamics nonrelativistically since as soon as the carriers leave the two 
dimensional world of the graphene sheet they attain mass (an adjustable parameter 
in the confinement model) and the hexagonal symmetry of the underlying lattice no 
longer applies indicating the nonrelativistic limit. This also indicates that the off- surface 
component of the wavefunction is a scalar which allows separation of variables. 

To achieve the above goal we will minimize the action of the Dirac equation ^ 



S= j dV ^[-ihvFY,l-Da + m{q^)vl + V^{q^) 



where the conjugated two component spinor wavefunction is \E' = '^^'-y^ . First, let us 
make few transformations. 

The volume element in the above action is defined with respect to a moving frame 
associated to the surface 

X{qi, q2, qs) = x{qi, ga) + QsNiqi, ^2), (9) 

— * 

where N(qi,q2) is a vector normal to the surface. This coordinate system facilitates 
the probing of the three dimensional space X{qi,q2,q3) from the viewpoint of a two 
dimensional surface x{qi,q2). For the volume element we obtain 

dV = fdSdq^, (10) 

where 

/ = 1 - 2Mg3 + Kql (11) 

and dS = ^/gdqidq2. Hereafter gij^s are the metric components and g is the determinant 
of the metric associated with the surface x(gi, ^2); M is the Mean and K the Gaussian 
curvature of this surface, respectively. 

In this construction the spinor wavefunction is normalized according to 

fdSdqs^-f^-^ = 1 (12) 

and in order to facilitate the separation of the dynamics into two parts, surface and 
orthogonal to it, we rewrite the action with the substitution 



Electronic properties of corrugated graphene: the Heisenherg principle 



6 



The normalization condition for the new spinor wavefunction ^jj is the following 

dSdq^i^j^ = 1 (14) 



and the action is 

^/fdSdqsip -ihvF^l°'Da + m{q3)vl + V±{q3) 

Now we pull the factor f'^^"^ through the gauge covariant derivative. For the 
corresponding equation of motion we obtain 

{-ihvFr"d3 + m{q3)vl + V±{qs)} ip 

-^n^pV7"f/'^a^ = 0, (15) 

a=0 V/ 

where is the remaining part of the gauge covariant derivative ([8| after the partial 



derivative ds has been removed. Here using the definition of / (11), 

f/^if-y^) = f/'M. (16) 

Now the tetrad fields, in the definition of the curved gamma matrices, constructed 
from the metric tensor ([9]) = e^rfg^ can be read off 

Wq = dt, (17) 

Wi = >/G^dqi, (18) 

W2 = \/Ch2dq2, (19) 

W3 = dqs, (20) 

where Gij are the corresponding metric components of ([9|. Note is constant which 
means that 7^ is constant. 

We are now ready to take into account the effect of the constraining potential 
Vx((?3) since when it applies the wavefunction is compressed between two steep potential 
barriers on both sides of the surface. The value of the wavefunction will be significantly 
different from zero only for very small range of values around ^3 = 0. Therefore, we take 
the limit gs — )■ in all coefficients in the above equation except the terms containing 
m(g3) and V±{q3). In this limit the tetrad constructed from the full metric tensor 
^^(Q'I) '?2, 93) — ^ eJJ(gi,g2) and the curved gamma matrices are defined with respect to 
the surface metric Gij — )■ gij 

7"(gi,g2)=e^(gi,g2)7^ (21) 
where the tetrad field Wa = ^adqii is given with respect to the surface metric gij 

Wo = dt, (22) 
^^1 = Vgudqi, (23) 
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W2 = y/g22dq2, (24) 
W3 = dq3. (25) 

In the limit we also have / — )■ 1 and T>3 — )■ —i-^A^. As a result we obtain the 
following equation 

{-ihvFi^d^ + m(g3)4 + ■^^(gs)} - l^^^^A^t/j 

c 

2 

- ihvF i^Daip - ihvF'y^Milj = 0. (26) 

a=0 

In Coulomb gauge ^3 = 0. The wavefunction is separable ip = XnXt into normal to the 
surface xn and tangential xt part. Since is constant, so is 7'^ and the equation is 
separable 

2 

ihvF-f^dtXT = -ihVF YDaXT ' thVFl'^MxT, (27) 

a=l 

ihvFl^dtXN = -ihvFj^ds + m{q3)v% + V^{q3)xN- (28) 

Here Da are gauge covariant derivatives defined with respect to the two dimensional 
metric of the surface x{qi,q2). The second equation in the case of graphene should be 
solved in the nonrelativistic limit since it describes the off surface dynamics. 

It turns out that the nonrelativistic system tends to extract more information from 
the manifold than the relativistic one since the relativistic fermion system only probes 
the surrounding space to first order in the derivatives of the wavefunction, whereas 
the Schrodinger fermions or bosons probe it to second order. More specifically, the 
Gauss curvature K appears only non-relativistically and is related through the Gauss- 
Bonet theorem to the topology class of the manifold. In the nonrelativistic limit this 
information is returned as higher derivatives of the wavefunction. As a result the Gauss 
curvature does not appear explicitly in the confining potential for the relativistic case. 



3. Graphene with a bump 

Now, let us explore a particular example, namely a bump on the surface of graphene. 
The bump has a Gaussian shape and is rotationally invariant. The reason we explore 
this geometry is the fact that it has been discussed previously [2] and we can compare 
our results which depend on the confining procedure with the results obtained by 
postulating the two dimensional Dirac equation to which the carriers of graphene are 
subjected. It turns out that due to the confining procedure an extra term appears which 
is proportional to the Mean curvature of the surface. 
In polar coordinates the surface is given by 

z{r) = Ae"'^'/'' (29) 
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Figure 2. A bump on the graphene surface affects local density of states (LDOS) as 



inferred by ( 46 1 . The shaded region is approximately where the Mean curvature is at 



maximum and the LDOS should be affected the most. 



and is asymptotically flat, lim^-^oo z{r-) = (see Fig. [2]). We will use this in order to 
check if the equation we obtain has the correct limit. Let us denote 



dz 



^-^^^ dr^ 



af{r)dr'^ 



\dr J 

which leads to the following for the line element 

ds"^ = dr^ + r^de^ + dz^ = [1 + a/(r)] dr^ + r^dO, 
where 9 G [0, 27r] is the polar angle. The metric is 



l + a/(r) 
r2 



(30) 



(31) 



(32) 



Before we derive the Dirac equation for the bump we will focus on the flat case 
which will provide the asymptotic. The fiat case is obtained by taking a = 0. Then the 
metric reduces to 




which yields one possible choice for the moving frame g^y = e"^el6ab '■ 



(33) 




The affine connection F^^ for this metric has the nonvanishing elements 



0, 



re 



which play a role in computing the spin connection [S] fl^ 



(34) 



(35) 



(36) 
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The spin connection takes part in the two dimensional Dirac equation on the metric (7^^ 
as follows 



trid^ + ^^.)^ = 0, 

where 7 = {'j^jVpY); see the Appendix. 
Finally the fiat hamiltonian is 



(37) 



H 



flat 



dr 



+ 



2r 



dr 



+ 



2r 







(3J 



Now that we know the asymptotic limit we may derive the corresponding 
hamiltonian for the Gaussian bump keeping in mind that due to confinement we have 
an extra term in the hamiltonian —ihvp'y^M as indicated by equation (27). This extra 
term in the 2 representation we are working reduces to 

V = -thvpM (39) 

due to the properties of the 7 matrices P, dUl [H] . In the curved geometry 



+ 



a 



dr 



ry/1 + af 27^(1 + a/)3/2 
where the principal curvatures of the surface. 



For the curved metric (32) the moving frame can be given by 



(40) 




which together with the afiine connection 



a 



dr 



yields for the spin connection 



J- 



l + af 



re 



1 



(41) 

(42) 
(43) 

2 a/1 + a/ 

Finally for the Dirac equation describing the carriers in the curved graphene sheet we 
have 



0, 
1 



1 



2 Vl + a/ 



[ir{d^ + n^)-ihvFM]ij = Q, 



(45) 



where 7 = (7°, vp"^^ iVp'^'^) is defined through the moving frame (21 ). The corresponding 
hamiltonian is 



bump 



+ 



-icr^dt — ihvpM 




(46) 



dr 







where the magnetic field associated with the gauge field 

^ = 1 
r 2rVl + «/ 



An 



(47) 
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is substantial 

= --drirAe) (48) 
r 

and for the typical size of the bumps on graphene varies between 0.5 - 3 T. It is 
straightforward to check the asymptotic form a = of this hamiltonian and to see 
that indeed 

Hbump — ^ Hfiat (49) 

when the bump flattens. 



One of the most obvious consequences of this geometric potential (39 ) is that it shifts 
the Fermi energy. Thus controlling graphene's electronic properties can be achieved not 
only by doping and creating imperfections in the lattice, but also by bending [T2]. 

4. Graphene in the context of a wormhole geometry 

Since we have developed an effective treatment of the massless carriers in graphene for 
the low-energy electronic properties of curved sheets in the continuum limit, we proceed 
to apply it to the wormhole geometry [T3]. Our motivation stems from one of the biggest 
gaps in modern physics, namely the relation between Quantum Mechanics and General 
Theory of Relativity. Constrained quantum systems present us with the opportunity to 
explore the quantum mechanical consequences of curvature. 
The metric of the wormhole is given by 

ds^ = -cMt^ + dl^ + (6^ + /2)(^02 ^ gi^2 ^^^2^^ ^5Q) 

where t is the proper time of a static observer, / = ± ^/r^ + 6q is the proper radial 
distance at constant time, h = h{l) is the shape function of the wormhole [6(0) = 6o is 
the radius of the throat of the wormhole] and (6',0) are spherical polar coordinates. 

As is usually the case with quantum systems we are interested in eigenmodes and 
we fix the time t = const. We also consider the case (p = 7t/2 which represents an 
equatorial cross-section of a three-dimensional wormhole at constant time. The line 
element therefore becomes 

ds^ = df + (6^ + f)de^, (51) 

which is precisely equivalent to the line element of the catenoid[T^ 

2 

ds^ = -^^dr^ + rH9\ (52) 
r'^ + Oq 

Note that at any arbitrary cross-section of the three-dimensional wormhole say (p = (po, 
the line element is 



2 

ds^ = -^^dr^ + eh^de\ (53) 
r"' + Oq 

where = sin^ G [0, 1]. For the catenoid this will mean only rescaling the radius from 
r to er. The catenoid with the biggest radius corresponds to the equatorial cross-section, 
see Fig. |3| 
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Figure 3. Any two-dimensional section (catenoid) of a three dimensional wormhole 
geometry. A large graphene sheet folded like this may serve as a solid state 
experimental realization of transport of relativistic particles trough a wormhole. In 
the process of folding the hexagonal symmetry of the underlying lattice should not be 
disrupted. 



In cylindrical coordinates {z, r, 6) a two dimensional cross-section of a wormhole is 
given by 



The principal curvatures ki and k,2 of this surface are 
1 z 

Ki = — sech^ — , K2 = —Ki. (55) 

This means that any two dimensional cross-section of a wormhole is a minimal 
surface 

M = ^(/ti + /€2) = 0. (56) 

Since the Fermi velocity in graphene is 0.3% of the velocity of light, the mixing 
of space and time coordinates is negligible. Therefore, the two dimensional catenoid 
section of the wormhole geometry can be realized through folded graphene providing us 
with a genuine solid state experiment. 

The constrained fermion dynamics we derived for the purposes of obtaining the 
electronic properties of curved graphene can now be combined with the geometry of the 
catenoid as an experimental solid state proposition testing the transmission of material 
particles through a wormhole. Indeed graphene can be shaped as catenoid provided the 
sheet is big enough so that while bending it one does not destroy the hexagonal symmetry 
of the underlying lattice and lose the relativistic description [5]. The destruction of the 
hexagonal symmetry can be inferred by the splitting of the G peak in the Raman spectra 
of graphene [151 IE]- 

It appears that due to the vanishing of the mean curvature, no extra potential 
subjects the relativistic particles heading for a transition "on the other side" of the 
wormhole. The governing equation is the free Dirac equation in curvilinear coordinates. 
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We have demonstrated that due to the essentially three dimensional nature of the 
carriers of graphene, using a purely two dimensional relativistic description to derive 
their properties contradicts the Heisenberg principle. In order to obtain the correct 
relativistic equation we employed a gradual confining procedure which produced a 
geometric potential in the equation of motion. As a result, the possibility to open a 
band gap in a purely mechanical way, simply by bending the material, is demonstrated. 
The result presented here also connects to the properties of quantum systems in the 
curved space-time of general relativity. It is interesting to note that due to the coupling 
of mechanical and electronic degrees of freedom through the geometric gauge potential 
(which for nanoscale radii of curvature is substantial) 

V = -ihvpM, 

mechanical vibrations of graphene can induce resonance electronic transitions, since the 
band gap is proportional to the curvature created by the mechanical oscillation itself. 
This can result in a new class of MEMS and NEMS based on the geometric properties 
of graphene. 

Even more intriguing is the observation that when a graphene membrane in 
MEMS and NEMS is set to oscillate this would result in a time-dependent Mean 
curvature M ~ sin {27rfrest) term, which can mechanically drive transitions between and 
within electronic bands [T7]. Strain induced modification of the electronic properties or 
"straintronics" is clearly an important field of study in graphene. Moreover, the band 
gap of the material can be tuned in the terahertz region where powerful sources of 
coherent radiation are much needed in industrial applications involving security aspects 
and medical scans. 

Appendix 

In 2+1 flat dimensions the fundamental 2 representation of the 7 matrices is satisfied 
by 

7o = -i(^3, It = -CTi 
for i = 1,2. Here 

"^=(? J)' "^=(! 7)' "^=(J -i)- 

They obey the following relations 

tr(7'^7%^) = -2e''^', 
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75 = 7o7i72, 

[75,7a] = 

Acknowledgments 

This work was supported in part by the U.S. Department of Energy. V.A. also 
acknowledges support by ICTP, Trieste, Italy where this paper was completed. 

References 

[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 
81, 109 (2009). 

[2] F. de Juan, A. Cortijo, and M. A. H. Vozmcdiano, Phys. Rev. B 76, 165409 (2007); A. Cortijo 

and M. A. H. Vozmcdiano, Phys. Rev. B 79, 184205 (2009) and references therein. 
[3] H. Jensen and H. Koppe, Ann. Phys. 63, 586 (1971). 
[4] R.C.T. da Costa, Phys. Rev. A 23, 1982 (1981). 

[5] F. Guinea, M.I. Katsnclson, and A.K. Gcim, Nat. Phys. 6, 30 (2010); V.M. Pereira, A.H. Castro 

Neto, and N.M.R. Peres, Phys. Rev. B 80, 045401 (2009). 
[6] S.A. Pulling, Aspects of Quantum Field Theory in Curved Space-Time, Cambridge University 

Press, Cambridge (1989). 
[7] E.-A. Kim and A.H. Castro Neto, Eur. Phys. Lett. 84, 57007 (2008). 

[8] S. Weinberg, Gravitation and Cosmology: Principles and Applications to the Ceneral Theory of 
Relativity, (Wiley, New York, 1972). 

[9] M. Burgess and B. Jensen, Phys. Rev. A 48, 1861 (1993). 
[10] A. Pnueh, J. Phys. A: Math. Gen. 27, 1345 (1994). 
[11] E. Flekk0y and J.M. Leinaas, Int. J. Mod. Phys. A 6, 5327 (1992). 

[12] V. Atanasov and A. Saxena, Phys. Rev. B 81, 205409 (2010); Y.N. Joglekar and A. Saxena, Phys. 

Rev. B 80, 153405 (2009). 
[13] M. Morris and K. Thorn, Am. J. Phys. 56, 395 (1988). 
[14] R. Dandoloff, A. Saxena, and B. Jensen, Phys. Rev. A 81, 014102 (2010). 

[15] T. M. G. Mohiuddin et al, Phys. Rev. B 79, 205433 (2009); S. Okada and A. Oshiyama, Phys. Rev. 
Lett. 91. 216801 (2003); Y. Miyamoto, S. Saito, and D. Tomanek, Phys. Rev. B 65, 041402(R) 
(2001); S. Piscanee, M. Lazzeri, J. Robertson, A.C. Ferrari and F. Mauri, Phys. Rev. B 75, 
035427 (2007); T. M. G. Mohiuddin et al, Phys. Rev. B 79, 205433 (2009). 

[16] S. Gupta and A. Saxena, J. Raman Spcctrosc. 40, 1127 (2009). 

[17] The frequencies fres of some modes in particular geometries have already been determined: D. 
Garcia-Sanchez et al, Nano Lett. 8, 1399 (2008). 



